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Abstract — The computation overhead is of major concern when 
going for increased accuracy in online power system security 
assessment (OPSSA). This paper proposes a scalable solution 
technique based on distributed computing architecture to 
mitigate the problem. A variant of the master/slave pattern is 
used for deploying the cluster of workstations (COW), which 
act as the computational engine for the OPSSA. Owing to the 
inherent parallel structure in security analysis algorithm, to 
exploit the potential of distributed computing, domain 
decomposition is adopted instead of functional decomposition. 
The security assessment is performed utilizing the developed 
composite security index that can accurately differentiate the 
secure and non-secure cases and has been defined as a function 
of bus voltage and line flow limit violations. Validity of 
proposed architecture is demonstrated by the results obtained 
from an intensive experimentation using the benchmark IEEE 
57 bus test system. The proposed framework, which is scalable, 
can be further extended to intelligent monitoring and control 
of power system 

Index Terms — first term, second term, third term, fourth term, 
fifth term, sixth term 

I. Introduction 

With the initiation of the deregulated electricity market, 
the system operators are concerned with the special measures 
to protect the system against severe contingences and to 
increase the security margins. These actions are performed 
by them based on the results obtained by conducting power 
system security analysis [l].The calculations required for 
the power system security assessment are performed based 
on the (n-1) criterion that requires the analysis of system 
behavior and the verification of operational limits violations 
for each credible contingency. Traditionally these analyses 
are carried out off -line as it requires the solution of system 
state equations in both static and dynamic time frame. These 
off-line analyses referred to as worst case scenarios, give 
operational limits often that are too restrictive or, in the case 
when the real time conditions differ to the reference values, 
highly conservative [2] . Therefore, these analyses appear to 
be inadequate in the new competitive scenario where there is 
an uncertainty in predicting the future operating conditions. 
This trend has increased the need for fast and more accurate 
methods of security assessment [l].The utilities are forced 
to conduct the real time power system security assessment, 
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in which an on-time estimation of system state variables is 
conducted using distributed data measurements [3] and the 
security is assessed in real time for a large set of probable 
contingencies and transactions. The real-time analysis could 
lead to a credible improvement of the utilization of the available 
infrastructure at adequate reliability levels allowing system 
operators to obtain more realistic operational guidance in 
planning preventive and corrective actions aimed to mitigate 
the effect of critical contingencies [1-2]. 

Traditional sequential computation is inadequate for on- 
line power system security analysis as the entire computation 
can take, typically less than a few minutes for the information 
to be useful [4]. The application of artificial intelligent [5] and 
probabilistic [6] based methodologies have been attempted 
for obtaining fast but less accurate solution for security 
assessment. 

Considerable research efforts [7]-[9] have also been 
oriented to develop dedicated computer architectures based 
on supercomputers or network of workstations for the fast 
solution of power system state equations. This method is 
applied particularly to on-line power system security 
assessment, where it is necessary to predict the impact of 
credible contingencies and suggest suitable preventive or 
corrective control actions within a few minutes to mitigate 
the effects of critical events. In recent years parallel 
processing based on distributed systems seems to be a viable 
solution to speed up the simulations in order to obtain results 
in useful time. Security constrained optimal power flow 
solution in a distributed computing environment is proposed 
in [8]. In [9] the various functions of security analysis are 
mapped on to a network of workstations which work as a 
continuous flow of base case conditions. As supporting tools 
in developing this activity, the application of TCP/IP based 
communication services and web based control architectures 
have been recently published in [2]. 

This work mainly focuses on power system static security 
assessment, contingency screening and ranking. Contingency 
screening and ranking is conventionally performed by 
computing a scalar performance index (PI), derived from DC 
or fast decoupled load flow solution for each contingency 
[10]. These methods generally employ a quadratic function as 
the performance index. This makes the contingency ranking 
prone to masking problems, where a contingency with many 
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small limit violations is ranked equally well with the one in 
which there are only a few large limit violations. Also, the 
selection of weighting factors in the performance index is 
found to be a difficult task, as it should be chosen based on 
both the relative importance of buses and branches and the 
power system operating practice [10]. In addition, majority of 
the performance indices do not provide an exact 
differentiation between the secure and non-secure states. 
The conventionally used performance indices were seen to 
be calculated separately for line flows and bus voltages, as 
the overall performance index defined as the sum or weighted 
sum of the scalar performance indices for bus voltages and 
line flows could not provide accurate results. 

In Ref. [10], authors have proposed an accurate method 
of critical contingency screening and ranking based on 

composite security index PI C which is calculated using 

Newton Raphson load flow technique. The PI is defined 

based on both bus voltage and line flow limit violations and 
it has been demonstrated in [ 1 0] that it completely eliminates 
the masking problem. It also provides a proper definition of 
security in which the secure state is indicated by an index 
value of '0' , while a value greater than ' 1 ' indicates an insecure 
state. Index values lying between '0' and ' 1 ' indicate the alarm 
limit. In this method, the difficult task of selecting the weights 
is also completely avoided. 

In this paper, a distributed computing architecture for 
on-line power system static security assessment based on 
composite security index is proposed and a prototype is 
designed. A variant of the master/slave pattern with only 
those tools in the open source domain are used for deploying 
the computational engine. Owing to the inherent parallel 
structure in security analysis algorithm, and to exploit the 
potential of distributed computing, domain decomposition is 
adopted. Experimental investigations are carried out on IEEE 
57 bus demonstrate the effectiveness of the proposed 
solution. 

The outline of the paper is as follows. Formalization of 
static security assessment problem is given in section II. 
Development of composite security index is given in section 
III. In section IV, a frame work for performing security 
assessment under distributed computing environment is 
presented. Deployment of computational engine along with 
the definitions of standard performance measures used in 
distributed computing are presented in section V. Experimental 
results and discussions are presented in section VI. Finally 
conclusions are drawn in section VII followed by references. 

II. Power System Security Assessment 

Power system security assessment is associated with the 
steady state and dynamic response of the power system to 
various disturbances. This process can be divided in to three 
sequential activities: i. contingency screening and ranking, 
ii. static and dynamic contingency analysis and iii. preventive 
and corrective control. The security analysis is performed 
according to the (n-1) criterion that requires systems to be 
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operated so as to withstand all single contingencies [1]. In 
this work the first and second activities are mainly considered 
as they are known to be the bottleneck in the online 
computations. 

A. On line Static Security Assessment 

The calculations needed for the on-line static security 
assessment requires the steady state solution of the power 
system state equations in order to identify the voltages in all 
network nodes and the power flows in each line in real time. 
This real time power flow solution, updated every few 
minutes, is adopted as reference in the automatic assessment 
of the static security of the system. The limit violations in 
bus voltages and line flows identified by computing a scalar 
performance index each for bus voltages and line flows. Then 
the solution engine automatically studies hundreds of 
possible contingencies that would happen on the power 
system determining how well the system can withstand them 
[2]. The sequence of major steps for on-line power system 
static security assessment is as follows: 

i. Acquire field data. 

ii. A software routine that solves the static power flow 
problem is invoked. This is then adopted in contingencies 
analysis as base case study for N configuration. 

iii. Check, if the network technical limits are violated. If 
violated the system is not secure in N configuration. 

iv. For each contingency, generate an input file containing 
the network data modified by the effect of the considered 
contingency. 

v. This file is then used by dedicated software routines to 
solve the corresponding power flow problem. 

vi. Check for each contingency, if the network technical 
limits are violated. 

vii. Generate alarms in the presence of an expected system 
malfunctioning. 

In this work, the violations in network technical limits are 
identified by computing for base case as well as for each 

contingency, the composite security index PI c proposed by 

the authors in [10], which is defined as a function of both 
power flow and bus voltage limit violations. Development of 
composite security index is discussed in the following 
section. 

III. The Composite Security Index 

In this paper, the composite security index PI C developed 

by the authors in [10] is used for static security assessment. 
The composite security index has two components one for 
bus voltage and the other for line flow security check. Two 
types of limits were defined for bus voltages and line loadings, 
namely the security limit and the alarm limit. The security 
limit is the maximum limit specified for the bus voltages and 
line flows. The alarm limit provides an alarm zone adjacent to 
the security limit, which gives an indication of closeness to 
limit violations. The alarm zone also provides a flexible means 
of specifying the cut-off point for contingency selection based 
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upon numerically ranked security index [10]. It is also possible 
to treat the constraints on the bus voltage and the line flows 
as soft constraints, thereby the violation of these constraints, 
if not excessive, may be tolerated for short periods of time. 
The system is considered insecure if one or more bus voltages 
or line flows exceed their security limit. If one or more bus 
voltages or line flows exceed their alarm limit without 
exceeding their security limit, the system is considered to be 
in the alarm state. If none of the voltages or line flows violates 
an alarm limit, the system is considered secure. This is 
indicated by an index value of '0' . 

It is assumed that the desirable voltage at each bus is 

known and is represented as vf . The upper and lower alarm 

limits and security limits of bus voltages are represented 

as F" , F. , v" , v\ respectively. The normalized upper and 

lower voltage limit violations beyond the alarm limits are 
defined asin(l). 



V. - F u 
d u = - — —■ if V > F" 

v,i , r d ' J i • 



v; 



d u v4 =o-, if v^f; 



F l -V. 
d ■ =— — • 

yd 



; if V^F' 
< i= 0; if Vj > F t l 



(1) 



For each upper and lower limit of bus voltages, the 
normalization factor g v i is defined in (2) 



u _ j J ; 



yd 



(2) 



v." 



For power flows, the limit violation vectors d p and the 

normalization factor g are defined in similar way.Since only 

the maximum limits are required to be specified for the power 
flow through each line, two types of upper limits are specified 

for each line, say the alarm limit P and the security limit P . 

The security limit is the specified maximum limit of the power 
flow through the line. The normalized power flow limit 
violation vectors for each line j can be defined as in (3). 



d , = if I P. I> P F , 

PJ Base MVA 1 Fj 

d pJ =0; if \P \<P, 



(3) 



where 17*, I is the absolute value of the power flow through 
the line j. 

The normalization factor for each line j, is defined in (4) as 
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For an N-bus, M line system, there are (N+M) dimensional 
normalized limit violation vectors of both bus voltages and 
line flows. In multi-dimensional vector space these limit 
violation vectors form a hyper -box and approximating the 
hyper-box by a hyper-ellipse inscribed within, a scalar valued 

index named as composite security index PI C [10] can be 
formed. The is defined in (5) as; 
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From the definition of composite security index, the system 
is said to be in one of the three states as follows. 

• Secure state if the PI C = 

• Alarm state if < PI C < 1 

• Insecure state if PI C > 1 

The contingencies can be accordingly ranked in 
descending order of severity based on. It is also possible to 
provide precise information about the buses and/or the lines 
in which the limit violations occurred so that proper control 
actions can be taken, without doing a detailed contingency 
analysis [10]. 

IV. Distributed Computing Framework for Performing 
Power System Security Assessment 

In the new competitive environment, in order to ensure 
high security and reliability levels in delivering electrical 
energy, there is a need to develop an accurate and high- 
performance security assessment method. To run the on-line 
security assessment algorithm described in the previous 
section, a framework based on distributed computing is 
proposed in this section, utilizing a hierarchical variant of the 
master/slave pattern that exploits the hierarchical topology 
of interconnected computers, such as a network of clusters 
that can assure high reliability, flexibility, high degree of 
scalability, and fault tolerance [11]. 

To reduce the execution time needed for security analysis, 
a concurrent algorithm based on domain decomposition is 
adopted, instead of functional decomposition that does not 
guarantee good performance [12]. In this method the whole 
job is divided in to similar tasks, each one assigned to a 
different processor. Since each task coincides with the 
sequential execution of the analysis of a single contingency, 
only minor modifications are necessary to the sequential 
algorithm discussed in section II. The activity diagram in Fig. 
1 shows the details of the adopted approach. The base case 
and a number of contingencies are analyzed to evaluate the 
security of the electrical grid. For each case a "slave" task is 
created, which does the power flow solution for the electrical 

grid, calculates the composite security index PI C using 
(5).The critical contingencies are then ranked according to 
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Fig. 1. Activity diagram for the Parallelized Security Algorithm 

severity based on, and eventually give alarms. 

V. Deployment of Computational Engine 

The computational engine is deployed on an architecture 
consisting of Intel ® Xeon ® CPU's operating at 2.33 GHz, 
4096 MB of static RAM. Each of the Xeon CPU is a quad 
core, implying for every CPU incorporated into the architecture 
there are 4 effective processors. All the CPU's are connected 
by a fast Ethernet local area network and are running on 
CentOS Linux. Only those tools in the open source are used 
to deploy the computational engine. A cluster of workstations 
is formed using the Rocks cluster distribution [13]. Thecluster 
of workstations (COW) can be further extended to cloud 
computing, as it is also equipped with Eucalyptus to facilitate 
deploying them as part of into Amazon EC2 based clouds. 
Ganglia is a scalable distributed system monitor tool for high- 
performance computing systems such as clusters and grids. 
It allows the user to remotely view live or historical statistics 
for all machines that are being monitored. 

A. Programming Tools Used 

Scilab equipped with Parallel Virtual Machine (PVM) or 
traditional C equipped with Message Passing Interface (MPI) 
is used for programming on the cluster. Scilab is a numerical 
computational package developed by researchers from the 
EMRIA and the Ecole nationale des ponts et chaussees (ENPC) 
[14]. The PVM computing model enables a collection of 
heterogeneous computer systems to be viewed as a single 
parallel virtual machine [15]. PVM transparently handles all 
message routing, data conversion, and task scheduling across 
a network of incompatible computer architectures. The 
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application is programmed as a collection of cooperating 
tasks. Tasks access PVM resources through a library of 
standard interface routines. These routines allow the initiation 
and termination of tasks across the network as well as 
communication and synchronization between tasks. The PVM 
software provides a unified framework in which parallel 
algorithms can be executed in an efficient and straightforward 
manner using existing hardware. MPI is a message-passing 
library interface specification [16]. MPI addresses primarily 
the message-passing parallel programming model, in which 
data is moved from the address space of one process to that 
of another process through cooperative operation in each 
process. MPI is a specification, not an implementation; there 
are multiple implementations of MPI. This specification is for 
a library interface; MPI is not a language, and all MPI 
operations are expressed as functions, subroutines, or 
methods, according to the appropriate language bindings, 
which for C, C++, Fortran-77, and Fortran-95, are part of the 
MPI standard. 

B. Performance Measures in Distributed Computing 

Distributed computing is the execution of a computer 
program utilizing multiple computer processors concurrently 
instead of using one processor. In its simplest form, the most 
obvious benefit of using this architecture is the reduction in 
execution time of the program. To measure the performance 
of the proposed algorithm on distributed environment, 
standard definition of two types of performance parameters 
viz. speed up factor, and computational efficiency [12]. The 
speedup refers to how much a distributed algorithm is faster 
than a corresponding sequential algorithm and is defined as 



T 

s n = — 

P rp 

P 



(6) 



where p is the number of processors, T ] is the execution time 
of the sequential algorithm on one processor and T p is the 
execution time of the parallel algorithm with p processors. S p 
therefore describes the scalability of the system as the 
number of processors is increased. Ideal speed up is p when 
using p processors, i.e. when the computations can be divided 
in to equal duration processes with each process running on 
one processor, with no communication overhead. The 
efficiency is a performance metric that describes the fraction 
of the time that is being used by the processors for a given 
computation. It is defined as 



P T p 



(7) 



It is a value, typically between zero and one, estimating 
how well-utilized the processors are in solving the problem, 
compared to how much effort is wasted in communication 
and synchronization. 

VI. Experimental Results 
The effectiveness of the proposed computational engine for 
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on-line static security assessment is demonstrated through 
experimental investigations carried out on IEEE 57 bus 
standard test system. The experiments involved the simulation 
of credible contingencies such as line outages for each test 
system for base load condition. Security assessment has been 
carried out as in [10] by computing the composite security 
index for the contingencies considered. For programming on 
the cluster, Scilab equipped with Parallel Virtual Machine 
(PVM) or traditional C equipped with Message Passing 
Interface (MPI) is used. 

For different contingencies, the composite security 
indices are computed as per (5) and are then ranked in the 
order of their severity. According to this index, the insecure 
cases are easily identified as those with values greater than 
'1' and the secure cases are defined as a value of '0' and can 
be excluded from the contingency list. If the index value is 
between '0' and '1', the system is in the alarm state. 

Experimental investigations are conducted by varying the 
number of processors used for computational engine in order 
to evaluate the reduction in the computation time for the 
analysis and the performance measures are analyzed. The 
results obtained for the test system are presented in the 
following sub-section. 

A. IEEE 57 Bus Test System 

The proposed method is applied to IEEE 57-bus test 
system. The system consists of 57 buses, 72 transmission 
lines and 8 transformers. The system data and single line 
diagram for IEEE 57 bus test system has been obtained from 
[17]. In order to define the composite security index, ±7% 
and +10% of the desired bus voltage is taken as the alarm 
limit and security limit respectively, for each bus. For the line 
flows, 80% of the thermal limit is chosen as the alarm limit 
[10]. 

The contingencies for which a security breach is observed 
are tabulated in Table I. Column 1 represents different line 
outage cases. For example, L 52-53 represents an outage of 
line connected between bus numbers 52 and 53. In Column 2, 

the composite security index PI c computed for the 

corresponding contingency case is presented and the 
security status is shown in column 3. In this case line outages 

L 52-53 and L 14-15 are found insecure for which the PI c 
value is greater than ' 1' . 

Table I. Contingency Ranking for IEEE 57 Bus Test System 
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For the line outages L 25-30, L 1-17, and L 5-6 the system is 
found to be in the alarm state with indices between '0' and 
'1'. For all other contingencies the system is found secure 
with an index value of '0' . The remaining contingency cases 
which are actually secure are not shown in the Table. 

The total execution time or turnaround time taken for a 
single processor architecture and multi-processor 
architecture with number of processors (P) in increments of 
two is tabulated in columns 2 and 4 of Table II, for both 
programming tools. The computational times are arrived at 
by performing the computation several times and the average 
has been taken. Percentage reduction in execution time in 
both cases is also analyzed and tabulated in columns 3 and 5. 

Table II. Computational Time for IEEE 57 Bus System 
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Scikb, PVM 


CMPI 


Execution 
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% reduction 
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time 
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^.364 


~9.~ 


i a 


7.096 


84.1 


:."j6 


84.3 
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It has been shown that, for IEEE 57-bus test system, the 
turnaround time required for complete static security assess- 
ment with a simulation engine based on single processor 
architecture is 44.639 and 36.232 seconds respectively for 
Scilab-PVM, and C-MPI programming tools. If multi-proces- 
sor architecture is used for computation, the turnaround time 
gets considerably reduced. For example with 6 numbers of 
processors the total execution time obtained is 12.682 sec- 
onds with a reduction in computation time of 71.6 % with 
respect to the single processor architecture for Scilab-PVM 
programming tool. 

Similarly, with C-MPI prograrriming tool, the execution time 
obtained is 10.382 seconds with a reduction in computation 
time of 71.3%. The performance of Scilab-PVM and C-MPI are 
compared by computing the performance parameters viz. speed 
up and computational efficiency as per (6) and (7) respec- 
tively. The variation of speedup and efficiency with the num- 
ber of processors, are given in Fig. 2 and Fig. 3 respectively. 

The speedup factor S p is increasing consistently with the 
increase in the number of processors. This implies that the 
performance improvement is guaranteed and hence the 
proposed system is scalable. It can be observed from Fig. 3 
that, as the number processors increased above 6 the 
computational efficiency also get improved. 

Conclusion 

In addressing the requirement for faster and accurate 
methodologies for real time power system security analysis, 
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a distributed computing based architecture is proposed 
utilizing the tools in the open source domain and the 
composite security index [10] is used for high speed and 
accurate security analysis for on-line applications. The 
preliminary investigations reveals that high performance 
computing engine based on COW's is a viable scalable 
solution for performing accurate and fast security analysis. 
The efficiency of the engine increases with the complexity of 
the system. As number of processors has increased, the 
computational cost gets decreased. This is a first step towards 
realizing the full potential of distributed architecture. The 
COW, being equipped with Eucalyptus, can be integrated 
into private or public clouds making use of the computing 
resources. As the cloud computing has the potential to furnish 
on-demand a dynamically variable computational power 
without affecting the accuracy, the research in this direction 
is under progress. 
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